perm filename QFACE.SAI[GEO,BGB] blob
sn#083868 filedate 1974-01-26 generic text, type C, neo UTF8
COMMENT ⊗ VALID 00002 PAGES
C REC PAGE DESCRIPTION
C00001 00001
C00002 00002 BEGIN "QFACE"
C00005 ENDMK
C⊗;
BEGIN "QFACE"
REAL AA,BB,CC,KK,QQ;
REAL X0,Y0,Z0;
REAL X1,Y1,Z1;
REAL X2,Y2,Z2;
REAL X3,Y3,Z3;
X1←+1; Y1← 0; Z1←0;
X2←-1; Y2←+1; Z2←0;
X3←-1; Y3←-1; Z3←0;
KK ← X1*(Z2*Y3-Y2*Z3) + Y1*(X2*Z3-Z2*X3) + Z1*(Y2*X3-X2*Y3);
AA ← (Z1*(Y2-Y3) + Z2*(Y3-Y1) + Z3*(Y1-Y2));
BB ← (X1*(Z2-Z3) + X2*(Z3-Z1) + X3*(Z1-Z2));
CC ← (X1*(Y3-Y2) + X2*(Y1-Y3) + X3*(Y2-Y1));
QQ ← 1/SQRT(AA*AA+BB*BB+CC*CC);
AA ← AA*QQ;BB ← BB*QQ;CC ← CC*QQ;KK ← KK*QQ;
OUTSTR(CVG(AA)&9&CVG(BB)&9&CVG(CC)&9&CVG(KK)&13&10);
OUTSTR(CVG(CC*1 - KK)&13&10);
END "QFACE";